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ABSTRACT 

Commissioning of SCUBA-2 included a program of skydips and observations of calibration sources intended to 
be folded into regular observing as standard methods of source flux calibration and to monitor the atmospheric 
opacity and stability. During commissioning, it was found that these methods could also be utilised to characterise 
the fundamental instrument response to sky noise and astronomical signals. Novel techniques for analysing on- 
sky performance and atmospheric conditions are presented, along with results from the calibration observations 
and skydips. 
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1. INTRODUCTION 

SCUBA-2 is a multi-wavelength submillimeter bolometric array receiver which is currently being commissioned 
on the James Clerk Maxwell Telescope (JCMT) at Mauna Kea, Hawaii. For full details of the instrument design, 
performance and commissioning details see the papers by Holland et al. EI (#7741-4), and Bintley et alP (#7741-5) 
in these proceedings. 

Calibration of astronomical observations in the submillimeter requires an accurate, high-frequency monitoring 
of the atmospheric opacity along the line-of-sight of the telescope and a measure of the instrument responsivity 
and transmission efficiency at the time of the observation. The methods used to calibrate SCUBA-2 data were 
based primarily on the successful methods of calibration used for SCUBA, and details of these methods and 
their results can be seen in Archibald et al.P Jenness et alP and Stevens & RobsonP Commissioning and the 
shared-risk observing carried out in December 2009 to March 2010 included a program of skydips to monitor 
the atmospheric opacity and stability and to investigate the transmission properties of the SCUBA-2 450/im and 
850/im bandpass filters. Nightly observations of astronomical sources with known flux properties were used to 
determine the instrument performance and to provide absolute flux calibration of the instrument. This paper 
presents the preliminary results from these commissioning observations. 



2. EXTINCTION CORRECTION AT SUBMILLIMETER WAVELENGTHS 
2.1 Atmospheric attenuation 

The atmosphere severely limits ground-based observations at submillimeter wavelengths, providing only a few 
semi-transparent windows even at a high, dry site such at the summit of Mauna Kea, Hawaii. SCUBA-2, like 
its predecessor, SCUBA (Holland et alP), has been designed to take advantage of the 850/Ltm and 450/xm atmo- 
spheric windows. The SCUBA-2 bandpasses at both wavelengths are shown superimposed on the submillimeter 
atmospheric transmission windows in Figure [l] The opacity in the submillimeter is primarily caused by wa- 
ter vapour, although both oxygen and ozone contribute. The transmission is wavelength dependent, with the 
shorter wavelength windows having lower transmission and faster susceptibility to worsening conditions. The 
water vapour column, and therefore the opacity of the atmosphere, can change on short timescales, and must be 
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Figure 1. 450/im (blue) and 850/im (red) bandpass filters superimposed on the submillimeter atmospheric transmission 
curve for Mauna Kea assuming 1mm of precipitable water vapour. 

monitored at high frequency in order to correct for these effects. 

Opacity conditions for night-time observing at the JCMT are assessed via the adjacent CSO 225GHz ra- 
diometer. The radiometer measures the sky opacity using a multiple-point skydip once every fifteen minutes. 
The JCMT has a 183GHz water vapour monitor (WVM) installed in the receiver cabin which measures the 
precipitable water vapour (PWV) along the line of sight at one second intervals^ . The PWV measured by the 
WVM is scaled with an empirical correction to produce a CSO-equivalent 225GHz opacity. The WVM and CSO 
values show excellent correlation for the stable part of the night between 9pm and 3am, with deviations occurring 
in the early evening and at sunrise when conditions become more variable. The r wvm to t cso relation is shown 
to b<P : 

T wvm = 1.01 * T cso + 0.01 (1) 

Previous investigations of the submillimeter opacity effects on bolometric observations at the JCMT^ , ref- 
erence the CSO opacity. This work will use the WVM measurements, using the above relation, as the WVM 
provides the time resolution necessary to analyse the high-frequency opacity changes that can significantly affect 
SCUBA-2 observations. 

2.2 Skydips 

The standard method for estimating the zenith sky opacity at a given wavelength and azimuth is by performing 
a skydip. SCUBA measured the sky brightness temperature directly as a function of elevation, with absolute 
temperature calibration achieved by observing a hot and cold load of known temperature!^ The SCUBA data 
was then fitted to a model describing both a plane-parallel atmosphere and the optical systemP 

Jrneas = (1 - Vtel)Jtel + VtelJatm ~ bw f T\ t el J atm exp(-TaU * A) (2) 

where J mea s is the measured sky brightness, Jtei and J a tm are the brightness temperatures of the telescope and 
atmosphere, respectively and r\ t ei is the transmission of the telescope. The final variable, bwf, is the bandwidth 
factor of the filter being used, and is a function of water vapour content. The sky temperature is not measured 
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directly by SCUBA-2. However, the detectors measures the change in optical loading with high accuracy. By 
measuring the change in optical load when opening a 1 K cold shutter located in front of the detector arrays an 
accurate measure of the optical load from the optics and sky is obtained. We ht this power change, AP/ leQt to 
the model. A full calibration of the relation between sky temperature and £±P heat was not completed during the 
initial commissioning run, however the response between the two is a linear relation. The data can therefore be 
fitted using a function of the form: 



AP heat =Pq-Px exp{-(a * {r wvm - (i)) * A) (3) 

This simple model does not take variation in temperature of the atmosphere but is accurate enough for the 
current purpose. A skydip is obtained by measuring the change in optical loading at a number of elevations. An 
example of fitting the equation to such a skydip is shown in Figure [2] 
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Figure 2. 450/mi (left) and 850/im (right) discrete skydip examples taken simultaneously at a PWV ~ 1mm. Square 
points are the heater power change APheat, as described in section [2~T| at each airmass, and the dashed line is the model 
fit to the data. 



2.3 Tau relations 

Until sufficient calibration data could be collected on the sky with SCUBA-2, it was assumed that the SCUBA-2 
filters as shown in Figure [T] were similar enough in their transmission properties to the SCUBA filters that the 
tau relations derived for SCUBA would be adequate for preliminary extinction correction. The original SCUBA 
opacity correction terms, as determined by Archibald et alP were Tsso = 4.05(r UJUm - 0.001) and T450 = 26.2(r UJt , m 
- 0.014) . 

Fifty skydips were completed in the period between February 22nd and March 14th, 2010. The individual 
dips were fitted to Equation [3j and an example of an individual fit of model parameters are shown in Figure [2j 
The advantage of the SCUBA-2 detector array design means that each time the shutter is opened from the 
fixed, closed shutter heater value, a new data point for APheat is obtained. Thus, every observation, not only 
skydips, can be used to determine the functional form of the relation between the sky power and the atmospheric 
transmission. If the wavelength-dependent r is accounted for correctly, the power to transmission relation should 
be linear. With the original SCUBA opacity corrections, the power-transmission data showed distinct curvature, 
indicating that these relations did not correctly account for the opacity as seen by the individual arrays. When 
the new relations were fitted to the skydips and then used to calculate the atmospheric transmission, of the form 




Figure 3. The crosses denote the 450/^m (left) and 850/xm (right) APheat for all observations between the 22nd of 
February and 14th of March 2010, as a function of atmospheric transmission, overlaid with the model fits to the skydip 
data (squares). 

Transmission = exp(-rxA), where r is the opacity and A is the airmass, the data showed an excellent linear 
correlation at both wavelengths. 

Figure [3] displays the results of the fit to Equation [3] The crosses in both figures indicate all APh eat for all 
observations taken between the 22nd of February and the 14th of March, 2010. Overlaid as green squares are 
the fitted data points from the fifty skydips obtained in that period. The final fitted form of the tau relations 
for each wavelength are described as: t 850 = 5.36(r wvm - 0.006) and t 450 = 19.04(r wt)m - 0.018). 

2.4 Astronomical calibrators 

Nightly flux calibration is achieved by observation of astronomical sources with known flux properties. Mars and 
Uranus were the primary calibrators, and were used predominantly in preliminary commissioning. During the 
observing period analysed here, the bulk of calibration observations were undertaken on a series of secondary 
calibrators, which are listed with their SCUBA fluxes in Table [I] 

Table 1. Secondary calibrators used for flux calibration of SCUBA-2. The flux values are sourced from the references 
noted in the table. 



Source 


RA(J2000) 


DEC(J2000) 


850/im flux 


450/xm flux 


HL Tau 


04 31 38.4 


+18 13 59.0 


2.36 ± 0.2<P 


9.9 ± 2.fP 


CRL 618 


04 42 53.60 


+36 06 53.7 


4.7 ± 0.37 s1 


12.1 ± 2.P 


CRT 2688 


21 02 18.81 


+36 41 37.7 


6.39 ± 0.5P 


30.9 ± 3.P 


IRC + 10216 


09 47 57.38 


+13 16 43.7 


8.8 ± l.P 


17.5+ 4.P 


V883 Ori 


05 38 19 


-07 02 2.0 


1.34 ± 0.01® 


7.28 ± 0.07® 


Alpha Ori 


5 55 10.31 


07 24 25.4 


0.628 ± 0.008® 


1.39 ± 0.0#! 


TW Hydrae 


11 01 51.91 


-34 42 17.0 


1.37 ± 0.01® 


3.9 ± 0.7® 


Arp 220 


15 34 57.21 


+23 30 09.5 


0.668 ± 0.007® 


2.77 ± 0.06® 



When an observation of a source is run through the map-making software it has been flat-fielded and extinc- 
tion corrected, and the resulting signal is in picowatts (pW). To convert the measured signal of a source from 
pW to Janskys, a 'flux conversion factor' (FCF) needs to be applied to the map. These FCFs are determined by 
regular observation of planets, when available, and the sources listed in Table [T] It is ideal to observe calibrators 
at least a few times over each night, to account for any variation in instrument performance and dish temperature 
changes over the course of the evenings observing. 



The FCF can be calculated in two ways. Firstly from the 'peak' signal of a point source, such that FCF pea k — 
Speak/ Ppeak where S pea k is the flux in the beam and P pea k is the measured peak power (pW) from the map, 
giving units of Jy beam -1 . Secondly the integrated value over an aperture can be calculated in units of Jy 
arcsec -2 , where FCF int = S to t/[PintA]. S to t is the total flux of the calibrator, P int is the integrated sum in pW 
and A is the pixel area in arcsec 2 . The method using the integrated flux was found to produce more consistent 
results as the peak value is susceptible to telescope effects such as focus and pointing drifts. Therefore, a third 
method was used for SCUBA-2 where the FCFi nt was taken and modeled with a Gaussian beam with a FWHM 
equivalent to that of the JCMT beam at each wavelength. The resulting FCF, denoted FCFbeamequiv) calculates 
a peak FCF from the integrated value assuming that the point source is a perfect Gaussian. The equation for 
FCF fceQmegtli „ is written as: 



S tot x 1.133 x FWHMl am 

Pint X A 



(4) 



where the FWHM(, eam for the 450/mi is 7.25" and the 850/im is 14.5", respectively. Ideally, when the instrument 
is stable and the flat-fielding and extinction correction are correctly accounted for, the FCF should be a constant 
factor. However, when the secondary calibrator "FGFbeamequiv values were calculated for the observing period 
from February 22nd to the 14th of March, and were plotted as a function of transmission, a trend could be 
observed that was particularly obvious at 450/im. This is shown by the green squares in Figure |4j When the 
new tau conversion factors calculated here are reapplied to the FCF's, the trend as a function of transmission 
is removed at both wavelengths, indicating that the new opacity corrections are approximating the extinction 
correction more effectively. 
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Figure 4. The 450/im (left) and 850/im (right) flux conversion factors (FCFs), in Jy/beam/pW, for all secondary calibrators 
observed in the commissioning period, as a function of atmospheric transmission. In both cases, the green squares indicate 
the FCF calculated with the SCUBA tau relation, and the black circles plot the FCF as calculated with the tau relation 
fitted to the SCUBA-2 skydips. The linear fits show that in both cases, the new tau relation successfully removes a 
significant transmission dependent term in the FCF data. 



3. SUBMILLIMETER SEEING 



Submillimeter seeing is caused by fluctuations in the refractive index of the atmosphere, largely caused by the 
passage of water vapour passing through the beam. The variations caused by this inhomogeneity are of particular 
concern for interferometers such as the Submillimeter Array (SMA), adjacent to the JCMT. The Smithsonian 
Astrophysical Observatory (SAO) operates a phase monitor which tracks the submillimeter seeing by measuring 
the path length shifts induced by the water vapour fluctuations in the atmosphere. A description of the SAO 
phase monitor can be found in Masson (1994) .ESI The phase difference is measured each second and then one 
minute's worth of data are analysed for the rms scatter to provide the seeing at 12GHz over three baselines. 

It was predicted that seeing would produce pointing offsets, however no correlation was seen between the 
SAO seeing and pointing excursions in SCUBAP Beam-broadening can occur if atmospheric turbulence on scales 
smaller than the telescope is significant. Experiments by Church and HillsEQ suggested that this was unlikely to 
be a large effect at JCMT. Attempts have been made previously to find a correlation between the seeing at the 
JCMT and the SMA seeing observed by the phase monitorP However, the standard SCUBA integration times 
were too long to allow for any direct measurement of the seeing to be reliably obtained, as seeing fluctuations 
instead occur on the order of a second or less. SCUBA-2, with a sampling rate of 200Hz, has the capability to 
image at these timescales. A recent data set taken on a night of extremely low tau (r wvm sa 0.03) yet extremely 
unstable localised weather prompted a close examination of the seeing using SCUBA-2 data. 



3.1 Mars 

The Mars observation shown in Figure [5] was obtained on the evening of March 10th, 2010. Early evening faults 
included a high rate of sky flat-field failures and inconsistent pointing results. The tau, though low, was extremely 
unstable. A mode of the SCUBA-2 reduction software allows the reduction of small time-slices of sequential data 
to be made into individual maps. With a source as bright as Mars, these maps can be made of only a few 
tenths of a second worth of data, and still have sufficient signal-to-noise for analysis. The Mars observation was 
a minute long, and individual maps span 0.4 seconds of data each. Figure [5] is a series of snapshots of these data 
over a minute, progressing in time from top to bottom. 

As can be seen from the image, both pointing shifts and beam distortion occur, and both effects can be seen 
to correlate between the two wavelengths. To determine the magnitude of the distortion caused by the seeing, 
the individual maps were analysed using the KAPPAE2I beamfit program. Figure [6] shows the measured centroid 
positions of Mars in each 0.4 second time-slice for a night (February 24th) of stable conditions {r wvm = 0.07) 
and March 10th (r wvm = 0.03). The left hand plot is typical of Mars observations taken throughout this time 
period. 

Eight SCUBA-2 observations of Mars taken between the 22nd of February and the 10th of March were anal- 
ysed using the method described above. The black circles and line in Figure [7] shows the rms of the 450/im 
FWHM of Mars in arcseconds. The green squares denote the rms phase of a single baseline of the SMA phase 
monitor for the same nights. The correlation between the SMA seeing and that observed with SCUBA-2 needs 
a more detailed analysis. The simple comparison in Figure [7] uses an average value for the SMA rms, calculated 
for the hour bracketing the time when the SCUBA-2 observation was taken. The trends in the two datasets 
show good agreement as the seeing deteriorates on the 9th and 10th of March. 



4. CONCLUSIONS 

SCUBA-2 calibration results from the commissioning and shared-risk observing period between February 22nd 
and March 14th of 2010 are presented and analysed. The atmospheric opacity relations for the 850/um and 
450/im bandpass filters are shown to be T85o = 5.36(r t0l)m - 0.006) and T450 = 19.04(r TO1)TO - 0.018), respectively. 
These relations were determined by an empirical fit to the skydip data obtained during this period. The relations 




Figure 5. The 450^im (left) and 850/im (right) Mars image (0.4 second samples) over the course of a minute's observation. 
Time progresses from top to bottom. 
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Figure 6. The 450/im offset centroid positions in RA and DEC (arcseconds) of Mars in 0.4 second slices for one-minute 
observations on the 24th of February (left) and the 10th of March (right). 
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Figure 7. The rms (arcseconds) of the 450/im FWHM (Mars) centroids for a series of observations between the 22nd of 
February and the 10th of March (black circles) versus the rms of a single selected baseline of the SMA phase monitor for 
the same nights (green squares). 



differ from those previously determined for the SCUBA instrument, which were being used in initial reductions 
of the science data taken with SCUBA-2. The flux conversion factors (FCF's) determined from observation of 
calibration sources were shown to have a trend as a function of atmospheric transmission at both wavelengths. 
Application of the new opacity terms in the extinction correction model removed this transmission dependence 
in the FCFs, thus improving the accuracy of SCUBA-2 calibration. 

An observation of Mars on the 10th of March 2010 was shown to exhibit the effects of submillimeter seeing. 
SCUBA-2's high sampling rate allowed time-slices of 0.4 seconds in the data to be analysed, and these sequential 
images showed evidence of both beam distortion and pointing shifts which were shown to lie well outside of the 
normal range when compared to similar observations on previous nights. 
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